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ABSTRACT 

> : 

}_j , We present the analyses of the stellar contents associated with the ex- 



tended H II region Sh2-252 using deep optical UBVRI photometry, slit and 
slitless spectroscopy along with the near-infrared (NIR) data from 2MASS 
for an area ^ 1 degree x 1 degree. We have studied the sub-regions of Sh2- 
252 which includes four compact-H ii (CH ii ) regions, namely A, B, C and 
E and two clusters NGC 2175s and Teutsch 136 (Ten 136). Of the fifteen 
spectroscopically observed bright stars, eight have been identified as massive 
members of spectral class earlier than B3. From the spectro-photomctric anal- 
yses, we derived the average distance of the region as 2.4±0.2 kpc and the 
reddening E{B — V) of the massive members is found to vary between 0.35 to 
2.1 mag. We found that NGC 2175s and Ten 136, located towards the eastern 
edge of the complex are the sub-clusters of Sh2-252. The stellar surface den- 
sity distribution in ii'-band shows clustering associated with the regions A, 
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C, E, NGC 2175s and Tcu 136. Wc have also identified the candidate ionizing 
sources of the CH ii regions. 61 Ha emission sources are identified using sht- 
less spectroscopy. The distribution of the Ha emission sources and candidate 
young stellar objects (YSOs) with IR excess on the V/{V — /) colour magni- 
tude diagram (CMD) shows that a majority of them have approximate ages 
between 0.1 - 5 Myr and masses in the range of 0.3 - 2.5 M©. The optical 
CMDs of the candidate pre-main sequence (PMS) sources in the individual 
regions also show an age spread of 0.1 - 5 Myr for each of them. We calculated 
the K-h&nd luminosity functions (KLFs) for the sub-regions A, C, E, NGC 
2175s and Ten 136. Within errors, the KLFs for all the sub-regions are found 
to be similar and comparable to that of young clusters of age < 5 Myr. We 
also estimated the mass functions (MPs) of the PMS sample of the individual 
regions in the mass range of 0.3 — 2.5Mq. In general, the slopes of the MPs 
of all the sub-regions are found comparable to the Salpeter value. 

Key words: stars: formation — stars: luminosity function, mass function — 
stars: pre— main— sequence — H ii regions: individual: Sh2-252 

1 INTRODUCTION 

It is known that most of the stars in the galaxy are formed in clusters or in OB associations 
where massive members of the region significantly affect their environment. This is possible 
in the form of strong stellar winds, intense UV radiation and eventually supernova explosions. 
These processes can ultimately cause their natal molecular cloud to be destroyed, which, 
in turn, will put an end to further star formation. However, it has also been observed that 
the above processes can trigger the birth of new generations of stars (e.g., Elmegreen & 
Lada 1977) and is likely to influence the key properties of the star formation process such as 
the initial mass function (IMF), star formation efficiency and the evolution of protostellar 
disks around the young stars (e.g., Clarke 2007). Keeping in mind the above points, it is 
important to study the large scale properties of OB associations so that we can have better 
understanding of the various processes governing the star formation. 

As a continuation of our multi- wavelength analyses of star forming regions (SFRs; e.g., 
Pandey et al. 2008, Jose et al. 2008; 2011), an optically bright H ll region Sh2-252 (Sharpless 
1959; 02000 = 06'^09™39'^; ^2000 = +20°29'12"; 1=190.04; b=+0.48) has been studied in this 
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paper. It is located in the galactic anti-center direction and is a part of the Gemini OBI 
association. Radio and CO surveys have already been done for this region by many authors 
(Felh et al. 1977, Lada & Wooden 1979, Kompe et al. 1989; Szymczak et al. 2000 etc.). 
However, stellar contents of this region was poorly studied in optical and infrared wavebands. 
Also, authors have raised arguments regarding the presence of many ionizing sources in this 
complex (see Felli et al. 1977, Lada & Wooden 1979). However, the massive members of 
this region have not yet been characterized spectroscopically except three bright sources. In 
this paper, we use the deep optical data in UBVRI bands, slit and slitless spectroscopic 
observations along with the JHK data from 2MASS to identify and classify the massive 
members, to determine the reddening, distance of the region as well as the age and mass 
of candidate pre-main sequence (PMS) sources and finally to obtain the K-band luminosity 
function (KLF) and IMF. After a brief description of the previous studies on Sh2-252, section 
2 describes various data sets used for the present study. Analyses and results including 
fundamental parameters of Sh2-252, identification and classification of massive stars in the 
associated clusters, its various properties such as age, distance, reddening, KLF and IMF 
are discussed in section 3 and section 4 summarizes the results. 

1.1 Overview of previous studies on Sh2-252 

The low density (ug = 9 cm~^) and the relatively large size (D ~ 20 pc) of the extended 
region place Sh2-252 in the class of evolved H ll regions (Churchwell 1974). Main source of 
ionization of Sh2-252 is the central star HD 42088 with a spectral type 06.5 V (Conti & 
Alschular 1971; Walborn 1972), which is a member of the Gemini OBI association. In Fig. 
[H the DSS2-R band image of the region around Sh2-252 for an area of ~ 1.3 x 1.3 square 
degrees is given, which shows a bright nebula having a diameter ~ 30' with a bright rim 
(Felli et al. 1977) at the west. See the colour composite image given in Fig. [2] for a close up 
view. Two stellar clusters were identified in this region optically, one is the dispersed cluster 
(namely NGC 2175) centered on the nebulosity and the other of smaller diameter (namely 
NGC 2175s) is located at the north-east of the nebula. Pismis (1970) estimated the distances 
to these clusters as 1900 pc and 3500 pc, respectively and hence NGC 2175s was thought 
to be a background object. Subsequently, using eight colour photometry, Chavarrfa et al. 
(1987) found that these two clusters are located at a distance of 2.0 kpc and hence claimed 
that both are associated with Sh2-252. A new near-infrared (NIR) cluster Teutsch 136 (here 
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Figure 1. The DDS2-R band image of tiie region around Sli2-252 for an area of ~ 1.3 X 1.3 square degrees. The locations of 
the thermal radio sources of the region A, B, C, E and F identified by Felh et al. (1977) along with the small clusters NGC 
2175s and Teu 136 are marked in the figure. The ionizing source (HD 42088) of the H II region is marked using a star symbol. 
The area covered for the deep optical observations in V and / bands is shown using continuous lines (polygon) and the thick 
dot-dashod box represents the area covered with the Ha slitless spectroscopy survey. The abscissa and the ordinates are for 
the J2000 epoch. 

after Teu 136; see Fig. [1]) has been identified towards the east of NGC 2175s by Koposov 
et al. (2008). However, they did not make any parameter determination for this cluster. 
Hence, the association of this object with Sh2-252 needs to be checked. In the literature, 
the distance to Sh2-252 from spectrophotometric and kinematic studies vary between 1.4 
kpc to 2.6 kpc (e.g., Bonatto & Bica 2011, Pismis 1970, Reid et al. 2009 and Grasdalen & 
Carrasco 1975). Using the 5 GHz aperture synthesis observations, Felli et al. (1977) detected 
six extended radio sources towards Sh2-252 labeled as Sh2-252 A to F (see Fig. 1 in Lada 
& Wooden 1979). They have identified four of these sources (A, B, C and E) as compact 
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Figure 2. A close-up colour composite view of Sh2-252 reproduced by using the V, Ha and [S ll] bands {V: blue; Ha: green; 
[S ll] : red) for an area ~ 50 X 50 arc min^ obtained from the Kiso observations. The sources numbered are further classified 
by using spectroscopic observations (see section [3. Il l and their identification are also shown in the figure. The locations of the 
ionizing source HD 42088 and the bright rim are also marked. The lines in the image are artifacts. 

H II (CH II ) regions, probably having a local source of ionization in each of them with a 
spectral type later than HD 42088 (i.e., 06.5V). The region F coincides with the optical 
bright rim (see Fig. [H and see also Fig. |2] for a close up view) and no internal heating source 
is found within it (Lada & Wooden 1979). The molecular survey by Lada & Wooden (1979) 
shows that the most intense CO peak is located very close to Sh2-252A with water and 
methanol maser emissions within its proximity (Lada & Wooden 1979; Kompe et al. 1989; 
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Szymczak et al. 2000), which indicates the recent star formation activity towards this region. 
Most of the region around Sh2-252A seems to be optically obscured, however, there is an 
embedded cluster visible in NIR and recently, Tej et al. (2006) has studied this cluster in 
detail. 

2 OBSERVATIONS AND DATA REDUCTIONS 
2.1 Optical CCD imaging 

The CCD UBVRI observations of Sh2-252 were carried out using the 105-cm Schmidt 
Telescope of the Kiso observatory, Japan, on 2004 November 4. The 2K x 2K CCD with 
each pixel corresponding to 1.5 arcsec covers a field of ~ 50 x 50 arcmin^ in the sky. The 
log of observations is tabulated in Table [1] 

The CCD UBVRI photometry of the central region (aaooo = OG^Og^SQ"; (52ooo = +20°29'12") 
of Sh2-252 was carried out using the 104-cm Sampurnand Telescope (ST) of ARIES, Naini- 
tal, India on 2005 December 2. The 2K x 2K CCD with a plate scale of 0.37 arcsec pixel"^ 
covers a field of ~ 13 x 13 arcmin^ in the sky. To improve the signal to noise ratio, the ob- 
servations were carried out in a binning mode of 2 x 2 pixels. The standard field SA 98 from 
Landolt (1992) was observed on the same night to apply the atmospheric and instrument 
corrections to the target field. We then performed PSF photometry with the DAOPHOT 
package in IRAF, and derived the extinction coefficients and colour coefficients using the 
observations of standard field SA 98. The central region of Sh2-252 observed with ST was 
calibrated by applying these coefficients. The calibration uncertainties between the standard 
and transformed V magnitudes and U — B — V ,V — R,V — I colours were ^ 0.03 mag. 
The selected isolated stars of this central calibrated region were used to calibrate the wide 
field observations taken with Kiso Schmidt Telescope. 

We further carried out deep observations in V and / bands using ST for 20 overlapping 
sub-regions of Sh2-252, to cover an area ~ 1 degree x 1 degree around it. The observations 
were conducted from 2006 October to 2008 December in good photometric conditions and 
the log of observations is tabulated in Table [H Each region was exposed for 60 min in V and 
30 min in /-band. The combined area of these 20 overlapping regions is shown in Fig. [1] using 
solid lines. The secondary standards from the Kiso observations were used to calibrate these 
individual regions. The final catalog is made from these three sets of observations, where, 
the magnitudes of the bright {V < 15 mag) and faint {V > 15 mag) objects were taken from 
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Table 1. Log of observations 



'-''{2000) '^(2000) Date of Filter & Exp. time (sec) X no. of frames 

(h:m:s) (d:m:s) observation 
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06:09:32 


+ 20 


28 


10 


2004 


11 


04 


U: 


60x2; B 




20x2; 


V: 10x2; R: 10 


































06:09:39 


+20 


29 


12 


2005 


12 


02 


U: 


360x2 


B 


: 240x2; V: 


120x2; R: 


06:10:09 


+20 


34 


57 


2006 


10 


19 


V: 


600x6 


I 


300x6 








06:09:14 


+20 


35 


60 


2006 


10 


23 


V: 


600x6 


I 


300 x( 










06:10:09 


+20 


25 


28 


2006 


10 


23 


V: 


600x6 


I 


300 x( 










06:09:18 


+20 


24 


52 


2006 


12 


06 
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600x6 


I 


300 xt 










06:08:06 
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20 


20 


2007 


01 


16 
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I 
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20 


20 
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01 


15 
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600x6 


I 


300 xt 










06:10:56 


+20 


35 


29 


2007 


01 


15 


V: 


600x6 


I 


300x6 








06:09:11 


+20 


25 


22 


2007 


11 


06 


V: 


600x6 


I 


300x6 








06:10:54 


+20 


36 


IS 


2007 


11 


06 


V: 


600x6 


I 


300 xt 










06:08:51 


+20 


12 


36 


2008 


01 


08 


V: 


600x6 


I 


300x6 








06:11:56 


+20 


38 


28 


2008 


01 


08 


V: 


600x6 


I 


300 xt 










06:10:25 


+20 


13 


46 


2008 


04 


15 


V: 


600x6 


I 


300 xt 










06:11:31 


+20 


34 


35 


2008 


04 


15 
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600x6 
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06:08:31 
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39 


24 
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26 
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I 
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46 


54 
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12 


19 
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I 
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450x3 




Hq 


Br 


90x2 


06:08:05 


+20 


40 


56 


2007 


01 


25 


Gr5/Ha-Br 




450x3 




Hq 


Br 


90x2 


06:08:43 


+20 


16 


58 


2007 


01 


26 


Gr5/Ha-Br 




450x3 




Hq 


Br 


90x2 


06:08:05 


+20 


25 


01 


2007 


01 


26 


Gr5/Ha-Br 




450x3 




Hq 


Br 


90x2 


06:08:05 


+20 


33 


01 


2007 


01 


26 


Gr5/Ha-Br 
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"'" 105-cm Schmidt Telescope, Kiso, Japan 

^ 104-cm Sampurnanand Telescope, ARIES, Naini Tal 

^ 2-m Himalayan Chandra Telescope, lAO, Hanle 



Kiso observations and ST observations, respectively. To ensure good photometric quality, 
we selected only those sources having uncertainty <0.2 mag in V^-band to make the final 
catalog. Thus we obtained photometry of 8791 sources detected at least in V and /-bands 
with a limiting magnitude of V ~23 mag. 

In order to check the photometry accuracy, we compared our photometry with the UBV 
photometries available in the literature by Grasdalen & Carrasco (1975) and Haikala (1994). 
The mean difference between our photometry and theirs were ^ 0.03 mag inV^B — V and 
U — B, which shows that our photometry is in agreement with the previous studies. 
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2.2 Ha slitless grism spectroscopy 

An Ha emission line survey of Sh2-252 was conducted between 2005 January and 2009 
November. The observations were done in the shtless mode with a grism as the dispersing 
element using the HFOSC of Himalayan Chandra Telescope (HCT) of lAO, Hanle, India. 
The central 2K x 2K pixels of the 2K x 4K CCD were used for data acquisition. The CCD 
with an image scale of 0.296 arcsec pixel"^ covers an area of ~ 10 x 10 arcmin^ in the sky. 
A combination of Ha broad-band filter (iJa-Br; 6100 - 6740 A) and Grism 5 were used. The 
resolution of Grism 5 is 870. We observed 25 overlapping sub-regions around Sh2-252 and 
the area covered for this observations is shown in Fig. [T] using a dot-dashed box. Multiple 
frames were taken to ensure the presence of Ha emitting sources. The log of the observations 
are given in Table [H Ha emission line stars with an enhancement over the continuum at the 
Ha line were visually identified. We obtained 61 Ha emitting sources within our surveyed 
area. The detection limit of the H-alpha survey is about 3 A in terms of equivalent width, or 
V ~ 22 in terms of magnitude. In Table [2], we have provided the J2000 coordinates, V^-band 
magnitudes and V — I colours for these emission line sources, if available. 

2.3 Slit spectroscopy 

Low resolution optical spectroscopic observations of fifteen optically bright sources in Sh2- 
252 were made using the HFOSC of HCT. These sources are the probable massive members 
of Sh2-252 and some of them are the candidate ionizing sources of the sub-H ll regions (see 
section [3.4. ip . The locations of these sources are marked in Fig. |2l The observations were 
taken on 16 and 17 November 2009. The spectra in the wavelength range 3800-6840 A with 
a dispersion of 1.45 A pixel"^ were obtained using the low resolution Grism 7 with a slit 
of 2 arcsec width and exposure time of 900s. One-dimensional spectra were extracted from 
the bias-subtracted and flat-field corrected images by using the optimal extraction method 
in IRAF. Wavelength calibration of the spectra was done by using an FeAr lamp source. 
Spectro photometric standard star (G191B2B) was observed on these two nights to give the 
flux calibration to the target spectra. 

2.4 Near-infrared data from 2MASS 

NIR JHK data for point sources within a radius of 30 arcmin around Sh2-252 have been 
obtained from 2MASS point source catalog (PSC) (Cutri et al. 2003). To improve photo- 
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Table 2. photometric data of Ha emission line sources in Sh2-252 



'^(2000) 
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err 


(V - /) 




06:07: 


;48 


.61 


+20; 


:39; 


:16. 
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- 
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06:07: 


:59, 


,15 


+20; 


:29; 


:19. 


.3 
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0. 


.053 
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.970 


0. 
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06:08 


:06 


,25 


+20; 


:33; 


:47. 


.0 


18, 
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1 
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06:08 


:30, 


,31 


+20; 
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" 




" 




06:08 
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0. 
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0. 
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06:08 


:33 


,26 
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.6 
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06:08 


:35 


.19 
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:22 
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- 
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:46 
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:39 
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06:08 
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0. 
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.120 





.012 


06:08 
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.3 
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.007 
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06:08 
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,77 


+20 


:36 
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metric accuracy, we used photometric quality flag (ph_qual = AAA) which gives a S/N ^ 
10 and a photometric uncertainty < 0.10 mag. This selection criterion ensures best quality 
detection in terms of photometry and astrometry (cf. Lee et al. 2005). The JHK data were 
transformed from the 2MASS system to the California Institute of Technology (CIT) system 
by using the relations given by Carpenter (2001). 
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2.5 Completeness of the data 

To study the luminosity function (LF) / mass function (MF), it is necessary to take into 
account the incompleteness of the data that could occur due to various factors (e.g., crowding 
of the stars). We used ADDSTAR routine of DAOPHOT to determine the completeness 
factor (CF). The procedure has been outlined in detail in our earlier works (see e.g., Pandey 
et al. 2001). Briefly, we randomly added artificial stars to both V and I images taken with ST 
in such a way that they have similar geometrical locations but differ in / brightness according 
to the mean {V — I) colour (~ 1.5 mag) of the data sample. The luminosity distribution of 
artificial stars was chosen in such a way that more number of stars were inserted towards 
the fainter magnitude bins. The frames were reduced using the same procedure used for the 
original frames. The ratio of the number of stars recovered to those added in each magnitude 
interval gives the CF as a function of magnitude. The minimum value of the CF of the pair 
(i.e., V- and /-bands) was used to correct the data incompleteness. It is to be noted that the 
spatially varying nebular background (see. Fig. [2]) and stellar crowding characteristics (see 
Fig. [9]) could affect the local completeness limit. In order to account for this, we estimated 
the CFs of two representative sub-regions Sh2-252C and NGC 2175s and the values obtained 
are given in Table El The CFs of both the regions seem to be agreeing well with each other 
except for the faint magnitude bin, where the variation seems to be significant. 

We also estimated the CFs of the 2MASS data using the i^-band images from archive 
around the regions Sh2-252C and NGC 2175s. Using the DAOPHOT / ALLSTAR package, 
we detected all the stars in the 2MASS PSC, with photometry accuracy better that 0.05 
mag. We then performed the completeness analysis using ADDSTAR routine as mentioned 
above. The CFs thus obtained for the two sub-regions are given in Table [3l The CFs of both 
the regions seem to be more or less similar in all the magnitude bins. 

3 RESULTS 

3.1 Spectral classification of the optically bright sources in Sh2-252 

The targets for the low resolution spectroscopy were selected on the basis of their brightness, 
their location within 2 arcmin radius of the sub-regions of Sh2-252 (see Fig. |2]) as well as their 
location on the optical colour-magnitude diagram (CMD; Fig. [7]). The coordinates of these 
sources are given in Table HI To determine the spectral type, we extracted the low- resolution, 
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Table 3. Completeness Factors (CFs) of optical and 2MASS /<-band data within the sub-regions Sh2-252C and NGC 2175s 





Optical 




2MASS i^T-band 




mag. range 


Sh2-252C 


NGC 2175s 


Sh2-252C 


NGC 2175s 


lU - 11 


1 nn 
l.UU 


1 nn 
l.UU 


1 nn 
l.UU 


1 nn 
l.UU 


11 - 12 


1.00 


1.00 


0.99 


0.99 


12 - 13 


1.00 


1.00 


0.97 


0.97 


13 - 14 


1.00 


1.00 


0.92 


0.95 


14 - 15 


1.00 


1.00 


0.85 


0.88 


15 - 16 


1.00 


1.00 


0.67 


0.70 


16 - 17 


1.00 


1.00 






17 - 18 


0.99 


0.99 






18 - 19 


0.96 


0.98 






19 - 20 


0.89 


0.92 






20 - 21 


0.80 


0.85 






21 - 22 


0.59 


0.71 







one dimensional spectrum of each source. The flux calibrated, normalized spectrum of these 
sources are shown in Fig. [3l 

In order to classify the stars earlier than B3, we used the criteria given by Walborn & 
Fitzpatrick (1990). For later type stars we used the criteria by Jacoby et al. (1984) and 
Torres-Dodgen & Weaver (1993). Spectra of O and B stars have the features of hydrogen, 
helium and other atomic lines (e.g., O ll, C ill. Si ill. Si IV and Mg ll etc.). Hydrogen and 
helium lines are usually seen in absorption for dwarfs whereas they may be in emission in 
super giants. In the case of early type stars, the ratio of He I 4471 /He ll 4542 is a primary 
indicator of the spectral type and the ratio is greater than 1 for spectral type later than 07. 
The line strength of He ll gets weaker for late O type stars and He ll 4686 is last seen in 
BO. 5 type stars (Walborn & Fitzpatrick 1990). If the spectrum displays He ll line at 4200 
A, along with the O ll/C ill blend at 4650 A, it indicates that the spectral type of the star 
must be earlier than Bl. The absence of He ll 4200, He ll 4686 and Mg ll 4481 and the 
weak features of silicon lines along with the weak O ll/C ill blends at 4070 and 4650 are 
the supporting criteria in the B1-B2 range. For type B2, He I is in its maximum and for 
later types Si ll 4128-4130 and Mg ll 4481 appears stronger (Walborn & Fitzpatrick 1990). 
The presence of the He I lines in absorption constrains the spectral types to be earlier than 
B5-B7. Presence of the spectral lines Mg II 4481 and Si ill 4552 is an indication of evolved 
early B type stars. The spectral lines such as Na I 5893, He I, Ha, Ca I 6122, 6162, Fe ll 
6456 are used to classify the late type sources (Torres-Dodgen & Weaver 1993). Finally, 
the spectral type was assigned to each source by visual comparison to the standard library 
spectra given in the literature (Jacoby et al. 1984) and our classification is given in Table HI 
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However, because of low resolution observations an uncertainty of ± 0.5 to 1 in the sub-class 
identification is expected. 

Eight out of the fifteen stars observed are found to have spectral types earlier than B3V. 
The sub-H ll regions of Sh2-252 (i.e., A, B, C and E) are found to have at least one star of 
spectral type early than B3V (i.e., star 7^1, #3, #5 and #11). It is quite interesting to note 
that the remaining four early type sources (i.e., star #6, ^^7, 7^8 and #9) belong to the small 
cluster NGC 2175s. The star ^2 shows Ha 6563 line in strong emission (not shown in Fig. 
|3]). This star also shows strong NIR excess emission in the (J — H)/{H — K) colour-colour 
(C-C) diagram (see Fig. E]) most likely originating from the extensive, hot circumstellar 
material around the star. The emission line spectra as well as the NIR excess emission are 
the primary indicators of an young stellar object (YSO) and hence we tentatively put this 
object in the Herbig Be category. 

The spectral types of some of these sources have already been mentioned in the literature. 
On the basis of the 1280 MHz radio continuum flux, Tej et al. (2006) predicted the spectral 
type of the ionizing source of Sh2-252A as BO-BO. 5V. The star 7^3, located within 2 arcmin 
of the radio peak of Sh2-252A, is the brightest source in the optical and NIR wavelengths. 
This could be the ionizing source of Sh2-252A. We classify this source as an 09.5V type 
star, which is in agreement with the classification by Tej et al. (2006). The star 7^5, the 
probable exciting source of Sh2-252E, had been assigned a spectral type between 08 to 
B0.5V by Garnier & Lortet-Zuckermann (1971) and BIV by Grasdalen & Carrasco (1975). 
Our classification (BIV) is in good agrrement with that given by Grasdalen & Carrasco 
(1975). The star #11, the brightest source of Sh2-252B was classified as B2V by Grasdalen 
& Carrasco (1975). However, our classification (BOV) leads to a class slightly earlier than 
their study. The star #2 had been identified as a probable Herbig Be source of spectral 
type B9 by Chavarria et al. (1989), which seems to match with our classification. Thus 
within errors the present spectral classifications are in accordance with those reported in 
the literature. 

3.2 Reddening, membership and distance of the bright sources in Sh2-252 

The extinction in a star forming region arises due to two distinct sources: (1) the general 
ISM in the foreground of the region [E{B — V)min], and (2) the localized cloud associated 
with the region [E{B - V)iocai = E{B -V)- E{B - V)„nn], where E{B - V) is the observed 
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Figure 3. Flux calibrated normalized spectra for the optically bright sources in Sh2-252. The star IDs and spectral classes 
identified are given in the figure. Important spectral lines are marked. 



reddening of the star embedded in the local cloud. The presence of the molecular cloud 
(Lada & Wooden 1979) and the nebulosity associated with Sh2-252 as evident in Fig. [2] 
would give rise to variable extinction within this region. To estimate the extinction towards 
the region, we used the {U-B)I{B-V) C-C diagram. Fig.lshows the (U-B)/{B-V) C-C 
distribution of all the stars detected in our Kiso observations, where the spectroscopically 
observed sources are numbered in red. It is to be noted that some of the spectroscopically 
observed stars are not covered/detected in the U and B bands, hence they are not included 
in Fig. m The star #16, which is the ionizing source of Sh2-252 (i.e., HD 42088), was 
saturated in our observations. Hence the data for this star is taken from Pismis (1977). 
The zero-age-main-sequence (ZAMS) by Girardi et al. (2002) that is reddened by E{B — V) 
= 0.23 mag along the reddening vector having a normal slope of E{U — B)/E{B — V) = 
0.72 seems to match with the first sequence of distribution in Fig. HI However, majority of 
the spectroscopically identified stars earlier than spectral type B3 show E{B — V) in the 
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Figure 3 continued. 



Table 4. Details of spectroscopically identified stars 



ID 


"(2000) 


^(2000) 


V 


B -V 




Spectral 


Vo - Mv 


D 


Associated 




(h:m:s) 


(d:m:s) 


mag 


mag 


mag 


class 


mag 


pc 


with 


1 


06:09:19.51 


+20:38:20.61 


11.79 


0.32 


1.66 


B2.5V 


12.15 


2690 


C 


2 


06:09:21.18 


+20:40:23.01 


14.11 


0.90 




B9Vc 






C 


3 


06:08:32.10 


+20:39:19.09 


13.88 




6.5 


09.5V 


11.63 


2120 


A 


4 


06:08:35.31 


+20:39:24.49 


12.50 


0.30 


1.28 


B8V 


11.47 


1970 


A 


5 


06:09:52.63 


+20:30:16.60 


11.10 


0.36 


1.92 


BIV 


12.37 


2980 


E 


6 


06:10:52.84 


+20:36:42.00 


11.05 




1.93 


B1.5V 


11.94 


2440 


N2175S 


7 


06:10:52.83 


+20:36:34.00 


10.34 




1.67 


BIV 


11.87 


2370 


N2175S 


8 


06:10:54.10 


+20:36:44.77 


11.51 


0.49 


2.27 


B2V 


11.69 


2180 


N2175S 


9 


06:10:55.97 


+20:36:29.65 


10.92 


0.49 


2.40 


B0.5V 


12.12 


2660 


N2175S 


10 


06:10:52.02 


+20:38:00.50 


10.95 


0.32 




F6V 






N2175S 


11 


06:09:00.14 


+20:38:28.14 


10.79 


0.55 


2.62 


BOV 


12.17 


2710 


B 


12 


06:09:06.49 


+20:39:01.43 


12.60 


0.17 


0.88 


B8V 


11.97 


2480 


B 


13 


06:09:00.65 


+20:35:04.74 


12.32 


0.74 




KOV 






B 


14 


06:11:46.38 


+20:39:21.10 


13.74 




2.77 


B6V 


11.87 


2360 


Tou 136 


15 


06:11:45.10 


+20:39:20.42 


15.85 






KOV-III 






Teu 136 


16* 


06:09:39.60 


+20:29:15.43 


7.55 


0.06 


1.14 


06.5V 


11.74 


2230 





For 16* data is taken from Pismis (1977) 



range 0.35 - 0.75 mag. The stars #10, ^^13 and #15 are found to be late type stars and 
hence may not be associated with Sh2-252. Assuming that uncertainties in the {U — B) 
and {B — V) colours mainly originate from their calibration (ref. section [2TT|) . the typical 
error in the estimation of E{B — V) would be ~ 0.04 mag. A careful inspection of the C-C 
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Figure 4. {U — B) / [B — V) colour-colour diagram for the stars in Sh2-252 region. The continuous curves represent the ZAMS 
by Girardi et al. (2002) shifted along the reddening slope of 0.72 (shown as dashed line) for E{B — V) = 0.23, 0.35 and 0.75 
mag, respectively. The spectroscopically observed sources are shown in red with their IDs as given in Fig. Eland in Table H 

diagram indicates the presence of further reddened population at {B — V) ~ 0.65 - 0.90 
mag. Some of them may be the reddened members of the Sh2-252 region. The population 
having E{B — V) ^ 1.00 - 1.25 mag could be part of the background Norma-Cygnus arm 
as discussed by Carraro et al. (2005) and Pandey et al. (2006). 

We also estimated the individual reddening E{B — V) of spectroscopically classified stars. 
We used the relation E{B — V)= {B — V) — {B — V)o, where {B — V) is the observed colour 
and {B — V)q is the intrinsic colour taken from Schmidt-Kaler (1982). Some of the stars do 
not have the {B — V) data (i.e., stars #3, #6, ^^7 and #14) and for these stars Ay was 
calculated by using the relations E{J — H) = {J — H) — {J — H)q and Ay = E{J — H)/0.11 
(Cohen et al. 1981), where, (J — H) colours were taken from 2MASS and intrinsic (J — H)o 
colours were taken from Koornneef (1983). The Ay thus calculated are given in Table HI 
The range of reddening of the spectroscopically identified stars shown in Fig. H] agree well 
with that given in Table HI 
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Figure 5. (J — H)/{H — K) colour-colour distribution of the spectroscopically classified stars in Sh2-252. Their IDs are as 
given in Table|4] The locus for dwarfs (thin solid curve) and giants (thick dashed curve) are from Bessell & Brett (1988). The 
red dotted lines represent the intrinsic locus of CTTSs (Meyer et al. 1997). The dashed line in black represents the reddening 
vector drawn from the early dwarf locus (Cohen et al. 1981). 

In Fig. |5]we have shown the (J — H) / (H — K) C-C distribution of all the spectroscopically 
identified stars in Sh2-252. The sources falling towards the right of the reddening vector are 
considered to have NIR excess (see sect. 13.4.31 for details). The star ^2 falls in the location 
where Herbig Ae/Be stars are generally found (Hernandez et al. 2005) and this justifies our 
spectral classification of this source as a Herbig Be type (see sect. 13. ip . The figure also shows 
that star ^^3 is highly reddened when compared to other massive stars in the field. 

The estimated spectral classification, apparent magnitudes and extinction of the massive 
members of Sh2-252 (Table H]) enable us to derive their distances. Using the present spectral 
classification along with Ay and V values from Table H] and My values from the table of 
Schmidt-Kaler (1982), we estimated the intrinsic distance modulus (Vq — Mv) of each source 
and the same is given in Table |H The average value of the intrinsic distance modulus is 
found to be 11.91 ± 0.18 mag, which corresponds to a distance of 2400ti8o PC- We also 
calculated the intrinsic distance modulus {Kq — Mk) of each source based on the observed 
colours from 2MASS data and intrinsic colours and absolute magnitudes from Koornneef 
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Figure 6. Vo/{B — V)o CMD for the spectroscopically identified massive member stars in Sh2-252. The continuous curve 
represents the ZAMS locus by Girardi et al. (2002) corrected for the distance of 2.4 Icpc. The sources are numbered as in Table 
|4] The average error in colour is shown in the lower left side of the figure. 

(1983) and Schmidt-Kaler (1982), respectively. The average value of the intrinsic distance 
modulus using 2MASS data was found to be 11.80 ± 0.10 mag, which is comparable with 
that of the intrinsic distance modulus estimated from optical data within errors. Fig. [6] shows 
dereddened Vo/{B — V)q CMD for the probable main sequence (MS) members of Sh2-252 
identified spectroscopically. In this figure, we have also plotted the theoretical ZAMS locus 
for solar metalicity {Z = 0.02) by Girardi et al. (2002), adjusted for the distance of 2.4 kpc, 
which seems to be matching well with the distribution of the probable MS members. The 
present distance estimate is in good agreement with that of Georgelin & Georgelin (1970; 
2.4 kpc), Neckel et al. (1980; 2.3 kpc) and Chavarria et al. (1987; 2.3 kpc). At the same 
time, it is also in agreement with the distance estimates by Grasdalen & Carrasco (1975; 
2.6 kpc) and Reid et al. (2009; 2.1 ± 0.03 kpc) within errors. 

An important issue to be discussed here is the association of NGC 2175s and Teu 136 
with Sh2-252. NGC 2175s is a small, optically visible cluster located at the east of Sh2-252 
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(at ^2000 = 06^10™52^; ^2000 = +20°36'36"; see Figs. □ and [2]). The photoelectric study by 
Pismis (1970) put this cluster at a distance of 3500 pc. However, in the present spectroscopic 
observations, we identified four early B type members of this cluster (i.e., stars #6, #7, #8 
and #9, see Table H]) and the distance estimates of these four sources are in good agreement 
with that of the average distance of Sh2-252. Hence, our spectroscopic analysis supports 
the association of NGC 2175s with the Sh2-252 complex. Similarly, the cluster Teu 136 was 
identified as a NIR embedded cluster by Koposov et al. (2008) in a survey of the 2MASS 
catalog designed to identify new open clusters. Located at coordinates 02000 = 06'^11™'56''; 
^2000 = +20°40'14", ~ 15' NE of NGC 2175s, this cluster lies close to the eastern edge of the 
Sh2-252 complex. We obtained spectra of two stars located close to this cluster (i.e., stars 
^14 and #15) and the distance of star #14 (see Table H]) seems to match well with that of 
Sh2-252. Also, the spectral types of two bright stars at the center of this cluster, which are 
saturated in our photometric observations are given as B1.5 and B2.5, respectively in Reed 
(2003). The average distance modulus of these two stars comes out to be 11.74 mag which 
seems to match with the average distance modulus of the Sh2-252 complex. Therefore we 
presume that Teu 136 could be associated with the Sh2-252 complex, too. Hence, in this 
paper we report that the NIR cluster Teu 136 as a sub-cluster of the Sh2-252 complex. 

In Fig. [71 we show the V/ (V—I) CMD for all the stars detected in our optical photometry. 
Spectroscopically identified massive MS members are shown with large filled circles in red, 
where the brightest star is HD 42088. The dashed curves are the ZAMS by Girardi et al. 
(2002) shifted for the distance of 2.4 kpc and reddening E{B -V)=0.35 and E{B-V)=0.75, 
respectively. The average photometric uncertainty in V magnitudes and (V — I) colours 
estimated from the difference between the input and recovered magnitudes of the artificial 
stars (see Sect. 12.51) are shown in the left side of the figure. The CMD clearly shows a 
broad MS band down to V ~ 14 mag with a significant number of sources located within 
the reddening strips (i.e., E{B — V) = 0.35 to 0.75 mag). In the absence of spectroscopic 
information, it is difficult to assign membership for these bright sources. Since Sh2-252 is 
located close to the Galactic plane (1=190.04; b=+0.48) at a distance of 2.4 kpc, it is quite 
natural that the region could be significantly contaminated by Galactic field MS and giant 
populations. The extinction variation inside the Sh2-252 region may cause a spread in the 
MS both for the member sources as well as for the background sources. A significant number 
of sources can also be noticed towards the right side of the ZAMS, part of which could be 
PMS stars. 
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Figure 7. V/ {V — I) CMD for all the stars detected in our optical photometry. The dashed curves are the ZAMS from Girardi 
ot al. (2002) corrected for the distance 2.4 kpc and reddening E{B — V) = 0.35 and E{B — V)=0.T6, respectively. The sources 
marked with red filled circles are the spectroscopically identified MS members of Sh2-252. The average error in {V — I) colour 
as a function of magnitude is shown in the left side of the figure. 

3.3 Age and mass distribution of candidate YSOs in Sh2-252 

In the absence of proper motion data, the membership of Sh2-252 can be estabhshed by 
identifying stars with indicators of youth. YSOs emit excess radiation in the infrared in 
comparison to MS stars due to the thermal emission from their circumstellar material. So 
YSOs can be identified on the basis of their infrared excess emission. We have also carried 
out a YSO survey of the Sh2-252 region using the 2MASS and the Spitzer-YRAC data sets 
(3.6, 4.5, 5.8 and 8.0 /im) and the detailed analyses will be presented in a forthcoming paper. 
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A total of 96 Class I and 400 Class II YSOs are identified within this complex on the basis of 
IR colour excess criteria given by Gutermuth et al. (2009). The possible contamination such 
as PAH-emitting galaxies, broad-line AGNs etc. have been excluded using the criteria given 
by Gutermuth et al. (2009). Similarly, the presence of H« emission line is considered as a 
significant characteristic of a YSO with ongoing disk accretion process (e.g., Dahm 2005). 
We identified 61 Ha emission line stars in our slitless spectroscopy survey of the Sh2-252 
region (see section I2l2|) . Our deep photometry provided V and / counterparts for 15 Class 
I, 179 Class II and 54 Ha emission line sources. Of these 54 Ha sources 4 are Class I and 33 
are Class II YSOs. The distribution of these candidate YSOs i.e.. Class I, Class II and Ha 
emission line sources on the V/ {V — I) CMD is shown in Fig. [H]by using red squares, green 
triangles and blue triangles, respectively. The distance to Sh2-252 is adopted as 2.4 kpc as 
estimated on the basis of the spectroscopic analyses mentioned in section 13.21 In Fig. |8] the 
ZAMS (thick solid curve) by Girardi et al. (2002) and the PMS isochrones (dashed curves) 
by Siess et al. (2000) for age 0.1 and 5 Myr are also shown. The ZAMS and isochrones 
are shifted for the distance of 2.4 kpc and average reddening i.e., E{B — V) = 0.5 mag. 
The average photometric uncertainty in V magnitudes and {V — I) colours estimated from 
the completeness simulation (see Sect. 12. 5p close to region C are shown in the right side of 
the figure. The distribution of YSOs on the V/ {V — I) CMD is an ideal tool to estimate 
the approximate ages of YSOs. It is evident from this figure that a majority of the YSOs in 
Sh2-252 are distributed between the PMS isochrones of age 0.1 and 5 Myr, which shows that 
the region Sh2-252 is significantly comprised of young sources. However, the photometric 
uncertainty, differential reddening, binarity, different evolutionary stages of YSOs etc. can 
cause a spread in the CMD. The spectroscopic analysis of the member stars shows that the 
reddening Ay in the region varies between 0.9 to 6.5 mag (see Table H]). Also, the NIR colour- 
colour analysis of the candidate PMS sources shows that they are reddened up to 8 mag 
(see Sect. I3.4.3p . which suggest that the region is significantly affected by the differential 
reddening. Similarly, the photometric uncertainties in the case of YSOs (see Fig. |H]) may 
be relatively large due to crowding and variable background which may cause the scatter 
in the distribution of the Class I and Class II sources in the CMD. Hence the differential 
reddening and photometric uncertainty would have an impact on the observed age spread 
in the CMD. The detailed spatial distribution and evolutionary status of the YSOs in the 
Sh2-252 region will be analyzed in the forthcoming paper. In Fig. [8] we have also shown the 
PMS evolutionary tracks by Siess et al. (2000) (thin solid curves) for various mass bins (the 
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Figure 8. V/(V - I) CMD for the identified candidate YSOs in Sh2-252 within the area shown by solid lines in Fig. [T] The 
Class I and Class II sources are shown with red squares and green filled triangles, respectively and the blue open triangles are 
the Hex emission line sources. The thick solid curve is the locus of ZAMS from Girardi et al. (2002), dashed curves are the 
PMS isochrones of age 0.1 and 5 Myr, respectively and the thin solid curves are the evolutionary tracks for various mass bins 
from Siess et al. (2000). All the isochrones and tracks arc corrected for the distance and reddening. The arrow indicates the 
reddening vector for Ay = 2 mag. The average error in (V — /) colour as a function of magnitude is shown in the right side of 
the figure. 



value of mass for each track is given towards its right), which indicate that a majority of 
the YSOs have masses in the range between 0.3 Mq - 2.5 Mq. 

The ages of the young clusters/H ll regions can also be derived from the post-main- 
sequence age of their most massive members. The most massive member of Sh2-252 is HD 
42088, an 06.5 star which is still in the MS (Conti & Alschular 1971; Walborn 1972) and 
is considered to be the main ionizing source of Sh2-252. Also, none of the spectroscopically 
identified OB stars of Sh2-252 are found to be evolved (see section [3?T|) . Hence, the age of 
the H II region Sh2-252 should be younger or of the order of the MS life time of the 06.5V 
star, i.e., ~ 4.4 Myr (Meynet et al. 1994). Thus we can put an upper limit to the age of 
Sh2-252 as ~ 4 Myr. However, this has to be taken as an approximate estimate as the low 
mass members of the region could have been formed prior to HD 42088. 
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3.4 Stellar contents of the sub-regions in Sh2-252 

3.4-1 Ionizing sources of the CH II regions 

As discussed in section [LT| the radio observation of this region by Felli et al. (1977) identified 
six compact thermal radio sources, of which the regions A, B, C and E have been claimed 
as CH II regions having its own ionizing sources in each of them. In order to search for the 
ionizing sources of these regions, we obtained spectra of optically bright sources within 2' 
radius around each region. The spectral types of these bright sources have been discussed 
in section 13.11 and each region is found to have at least one source with spectral type earlier 
than B3V. Their spectral type can also be obtained by measuring the number of Lyman 
continuum photons (N/yc ) emitted per second by the star that would be responsible for the 
ionization of each region. From the radio continuum flux at 1415 MHz by Felli et al. (1977) 
together with our estimated distance of 2.4 kpc and assuming an electron temperature of 10"^ 
K, we calculated log N^^c using the relation given by Mart in- Hernandez et al. (2003). The log 
^lyc have been obtained as 46.73, 47.52, 46.57 and 46.55, respectively for regions A, B, C and 
E. From these log N;,yc values we derived the spectral type of the ionizing sources (Panagia 
1973) as BO-BO. 5, BO, BO. 5 and BO. 5, respectively, for regions A, B, C and E. These spectral 
types are consistent with the earliest spectral types derived by our optical observations for 
regions A, B and E (i.e., star #3, ^^11 and #5, see TableH]) within the errors. The 1280 MHz 
observations by Tej et al. (2006) also estimated the spectral type of the ionizing source of 
region A as BO-BO. 5V which is consistent with our estimate. Slight difference in the spectral 
type estimate is expected due to uncertainties in the effective temperature and ionizing flux 
of massive stars, for a given spectral type. Thus we conclude that the stars ^^3, #11 and #5 
are the ionizing sources of the CH ll regions A, B and E, respectively. However in the present 
observation, the earliest spectral type derived in region C is B2.5V (see Table Hj), which does 
not seem to be early enough to explain the radio continuum flux. The expected spectral type 
from the radio flux is much earlier than B2.5V. The study by Kompe et al. (1989) showed 
that the probable spectral type of the IRAS point source (IRAS 06063+2040) located at the 
center of the nebulosity in region C is BOV, which could be the candidate ionizing source 
of this CH II region. Since this source is deeply embedded in the cloud (i^=9.16. Ay ~ 8 
mag), we could not carry out optical spectroscopy of this source and future spectroscopy 
would help us to conflrm this result. In Fig. [9] we have shown the /T^-band images of the 
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Figure 9. 8' x8' ii'g-band images of the sub-regions within Sh2-252. The contours arc the stellar surface density distribution 
generated by using all the sources detected in X^-band. The ionizing sources of the CH II regions are also marked using white 
star symbol. 

sub regions A, B, C, E, NGC 2175s and Teu 136 for an area 8' x 8'. The probable ionizing 
sources of the regions A, B, C and E are marked with a white star symbol. 

3.4-2 Stellar surface density distribution in Sh2-252 

Though we have the optical and NIR observations for an area ~ 1 degree x 1 degree around 
the Sh2-252 region, however a majority of the members of Sh2-252 might be confined within 
its sub-regions (see section II. ip . In order to define the extent of each sub-region and to 
study the morphology of the region, we generated isodensity contours for stellar population 
detected in 2MASS i^-band. The two dimensional stellar surface density distribution of the 
Sh2-252 region manifests prominent stellar density enhancements at the locations of A, C, 
E, NGC 2175s and Teu 136. This indicates the presence of clustering associated with these 
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sub-regions. On the other hand, no clustering is apparent towards the sub-regions B and 
F. In Fig. M the stellar surface density contours above 3a of the background level are over 
plotted on the i^5-band images of regions A, B (no contours although), C, E, NGC 2175s and 
Teu 136 . We will discuss the properties of these five sub-regions i.e.. A, C, E, NGC 2175s 
and Teu 136 individually in ensuing sections. The extent of these sub-regions is estimated 
on the basis of the surface density distribution. The contour above the 3a background level 
around each region is considered as its extent. Thus we estimated the radius as 3'. 5, 3'. 5, 
3'.0, 3'.0, 3'.0, respectively, for the region A, C, E, NGC 2175s and Teu 136. We used these 
extents for further analysis of these sub-regions. 

3.4-3 NIR colour-colour diagrams 

NIR C-C diagrams are ideal tools to identify the candidate PMS sources having NIR excess 
(Hunter et al. 1995; Haisch et al. 2000; 2001; Sugitani et al. 2002; Devine et al. 2008; 
Chavarria et al. 2010). The (J — H)/{H — K) C-C diagrams for the sub-regions A, C, E, 
NGC 2175s and Teu 136 and a nearby control field are shown in Fig. [TOl We have considered 
an area of 3'.0 radius centered at 02000 = 06'^09"30^ ^2000 = +20°08' as the control field. 
This region is devoid of nebulosity in the i?-band (see Fig. [1]), hence can be considered as a 
field region. The thin and thick solid curves are the locations of unreddened MS and giant 
stars (Bessell & Brett 1988), respectively. The dotted lines represent the locus of intrinsic 
and reddened {Ay = 4.0 mag; 8.0 mag) classical T Tauri stars (CTTSs; Meyer et al. 1997). 
The two long parallel dashed lines are the reddening vectors for the early MS and giant 
type stars (drawn from the base and tip of the two branches). One more reddening vector is 
plotted from the tip of the unreddened CTTS locus. The crosses on the reddening vectors are 
separated by an Ay value of 5 mag. The extinction ratios, Aj/Ay = 0.265, A^/Ay = 0.155 
and Ax/Ay = 0.090, are adopted from Cohen et al. (1981). The magnitudes, colours of the 
stars and the curves are in the CIT system. 

Following Ojha et al. (2004a), we classified sources according to their locations in (J — 
H)/{H — K) C-C diagrams. The 'F' sources are those located between the reddening vectors 
projected from the intrinsic colours of MS and giant stars. These sources are reddened 
field stars (MS and giants) or Class Ill/Class II sources with little or no NIR excess (viz., 
weak-lined T Tauri sources (WTTSs) but some CTTSs may also be included). The sources 
located redward of region 'F' are considered to have NIR excess. Among these, the 'T' 
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Figure 10. {J — H)/{H — K) C-C diagrams of the sources detected in the J H K-hands within the sub-regions of Sh2-252 and a 
nearby control field. The regions are marked in the figure. The locus for dwarfs (thin solid curve) and giants (thick solid curve) 
are from Bessell &; Brett (1988). The red dotted lines represent the intrinsic and reddened {Ay = 4.0 mag; 8.0 mag) locus of 
CTTSs (Meyer et al. 1997). The dashed lines in black represent the reddening vectors (Cohen et al. 1981). The crosses on the 
dashed lines are separated by Ay = 5 mag. The blue vertical line divides the {H — K) colour at 0.35 mag (see the text). 



sources are located redward of 'F' but blueward of the reddening line projected from the red 
end of the CTTS locus. These sources are considered to be mostly CTTSs (Class II objects) 
with large NIR excesses (Lada & Adams 1992). There may be an overlap in NIR colours 
of Herbig Ae/Be stars and T Tauri stars in the 'T' region (Hillenbrand et al. 1992). The 
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'P' sources are those located in the region redward of region 'T' and are most hkely Class 
I objects (protostellar-like) showing large amount of NIR excess. Here it is worthwhile to 
mention that Robitaille et al. (2006) have shown that there is a significant overlap between 
protostellar-like objects and CTTSs in the C-C diagram. 

A comparison of the distribution of sources in the sub-regions and the control field in 
Fig. [To] suggests that there is an appreciable difference between them. A significant fraction 
of sources within the sub-regions are concentrated between the intrinsic and reddened CTTS 
locus (i.e., within the 'F' and 'T' regions), whereas a majority of the sources in the control 
field are mainly concentrated towards the left of the 'F' region (i.e., [H — K) < 0.35 mag). 
Average value of the uncertainty in {H — K) colour is ~ 0.05 mag. Taking this uncertainty 
into consideration and on the basis of the comparison of the C-C diagrams, we can safely 
assume that a majority of the sources in the sub-regions located between the intrinsic and 
reddened CTTS locus with [H — K) > 0.40 mag and lying in the 'F', 'T' and 'P' regions 
are most likely PMS members. Some of the above sources lying in the 'F' region could be 
the reddened field stars, however a majority of them are likely candidate WTTSs or CTTSs 
with little or no NIR excess. Fig. [TO] shows that the clusters NGC 2175s and Teu 136 do not 
have many PMS members with Ay > 4.0 mag, whereas a significant fraction of the PMS 
members of the sub-regions A, C and E are reddened up to 8.0 mag. This indicates that the 
clusters NGC 2175s and Teu 136, which are located towards the east of Sh2-252, are less 
reddened in comparison to other sub-regions of Sh2-252 located towards the west. 

3.4-4 Optical colour-magnitude diagrams 

The V/{V — I) CMDs for the sources within the sub-regions of Sh2-252 and the nearby 
control field are shown in Fig. [TT] The ZAMS (thick solid curve) by Girardi et al. (2002) 
and PMS isochrones (dashed curves) by Siess et al. (2000) for age 0.1 and 5 Myr are also 
shown. The ZAMS and isochrones are shifted for the distance of 2.4 kpc and reddening. The 
encircled sources are those located above the intrinsic CTTS locus in the (J — H) / {H — K) 
C-C diagram (Fig. [TO]) and with {H — K) > 0.4 mag. A majority of these sources are found 
to have ages < 5 Myr in all the sub-regions, suggesting that the sources with {H — K) > 
0.4 mag could be probable PMS members as discussed in section 13.4.31 A comparison of 
these CMDs with that of the control field also supports the notion that the sources lying 
between 0.1 and 5 Myr isochrones could be the PMS sources associated with the sub- regions. 
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Figure 11. V/{V — I) CMDs for the sources within the sub-regions of Sh2-252 and the nearby control field. The sub-regions are 
marked in the figure. Encircled are the sources with {H — K) > 0.4 mag and are located above the intrinsic CTTS locus of the 
(J - H)/{H - K) C-C diagram (Fig.[TOll. The thick solid curve is the ZAMS from Girardi ct al. (2002) and the dashed curves 
are the PMS isochrones of age 0.1 and 5 Myr, respectively, from Siess et al. (2000), corrected for the distance and reddening. 



Similarly, the CMDs of the sub-regions resemble the V/ {V — I) CMD for the candidate YSOs 
shown in Fig. [Hi It further supports that the majority of the sources identified on the basis 
(J — H)/{H — K) C-C diagrams (Fig. [TOl) in the sub-regions could be the PMS members 
and that they have an age spread of 5 Myr. A detailed analysis on the evolutionary status 
of the YSOs within the individual regions will be done in the forthcoming paper. 
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3.5 i^-band luminosity functions of sub-regions in Sh2-252 

The luminosity function in the K- band is frequently used in studies of young clusters and 
star forming regions as a diagnostic tool of the mass function and the star formation history 
of their stellar populations. Pioneering work on the interpretation of KLF was presented by 
Zinnecker et al. (1993). During the last decade several studies have been carried out with the 
aim of determining the KLF of young clusters (e.g., Muench et al. 2000; Lada & Lada 2003; 
Ojha et al. 2004b; Sanchawala et al. 2007; Pandey et al. 2008; Jose et al. 2008; 2011). We 
used the 2MASS K- band data to study the KLF of the Sh2-252 region. Because this region 
shows clusterings around the regions A, C, E, NGC 2175s and Teu 136 (see section [3.4.2p . 
we estimated the KLF for each region independently. The KLF is calculated within the area 
of each region as estimated on the basis of stellar surface density distribution discussed in 
section 13.4.21 

KLF is the number of stars as a function of K- band magnitude. In order to convert the 
observed KLF to the true KLF, it is necessary to account for the data incompleteness as 
well as the background and foreground source contamination. The regions C and NGC 2175s 
were corrected for their respective CFs given in Tableland the average value of the CFs for 
regions C and NGC 2175s was applied to regions A, E and Teu 136 to correct for their data 
incompleteness. The data from the control field region (see section I3.4.3P has been used to 
remove the field star contribution. Since this control field is located off the Sh2-252 region, 
its background population will be affected by a smaller interstellar reddening. Whereas, the 
background field stars of the sub-regions of Sh2-252 are seen through a larger reddening due 
to the matter associated with it. Hence the population in the background of sub-regions will 
be more reddened as compared to the control field. Therefore a direct subtraction of the 
observed control field from the sub-regions of Sh2-252 will yield incorrect LF. To account for 
the higher extinction towards the embedded sub-regions, the field star population towards 
the direction of the control field is predicted on the basis of the Besangon Galactic model of 
stellar population synthesis (Robin et al. 2003) by using a similar procedure as described by 
Ojha et al. (2004b). An advantage of using this model is that we can simulate foreground 
(d < 2.4 kpc) and the background (d > 2.4 kpc) field star populations separately. The 
foreground population was simulated by using the model with the extinction Ay = 1.1 
mag {E{B ~ V) = 0.35 mag; cf. section 13.21) and d < 2.4 kpc. As discussed in section 
I3.4.3[ a majority of the PMS members of regions A, C and E are reddened up to Av=8.0 
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mag, whereas in NGC 2175s and Teu 136 they are reddened up to y4v'=4.0 mag. Hence the 
background population in these regions are seen through clouds with extinction of 8.0 and 
4.0 magnitudes, respectively. Hence, we simulated two sets of background population (d > 
2.4 kpc) with Ay values of 4.0 and 8.0 mag, respectively. Then we determined the fractions 
of the contaminating stars (foreground + background) over the total model counts. The 
fraction was calculated separately for both sets of the simulated background population. 
This fractions were used to scale the nearby observed control field and subsequently the star 
counts of the modified control field were subtracted from the KLFs of the sub-regions of 
Sh2-252 to obtain their final corrected KLFs. 

The KLFs of young embedded clusters are known to follow power-law shapes (Lada et 
al. 1991; 1993), which is expressed as: 

dN(K) . r.aK 
dK 

where is the number of stars per 1 mag bin and a is the slope of the power law. 

We counted the number stars per unit magnitude interval after correcting for the field star 
contribution. We then fit a power law to the data. The resulting field star subtracted KLFs 
for the five sub-regions of Sh2-252 are shown in Fig. [121 The slopes of the KLFs of regions 
A, C, E, NGC 2175s and Teu 136 have been obtained as 0.31±0.03, 0.34±0.03, 0.25±0.05, 
0.20±0.05 and 0.36±0.02, respectively. Within errors, the KLFs for all the five regions seem 
to match with each other. 

KLFs of different ages are known to have different peak magnitudes and slopes (Muench 
et al. 2000) and hence the KLF slope could be an age indicator of young clusters. For clusters 
up to 10 Myr old, the KLF slope gets steeper as the cluster gets older (Ali & Depoy 1995; 
Lada & Lada 1995). However, there exists no precise age - KLF relationship in the literature 
due to huge uncertainties in their correlation (Devine et al. 2008). There are many studies on 
KLF of young clusters. The studies by Blum et al. (2000); Figueredo et al. (2002); Leistra 
et al. (2005; 2006) and Devine at al. (2008) indicate that the KLF slope varies from 0.2 
- 0.4 for clusters younger than 5 Myr, which is in agreement with the KLFs obtained for 
the sub-regions of Sh2-252. The KLFs of the sub-regions of Sh2-252 are worth comparing 
with the recent studies of star forming regions viz; NGC 1893 (Sharma et al. 2007), Be 59 
(Pandey et al. 2008), Stock 8 and NGC 1624 (Jose et al. 2008; 2011), since all the KLFs are 
obtained by using a similar technique. The slopes of the KLFs obtained for the sub-regions 
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Figure 12. KLFs derived for the sub-regions in Sh2-252 after subtracting the field star contamination (see the text). The error 
bars represent the ±\/]V errors. The linear fits are represented by the straight lines and the slopes obtained are given in each 
figure. 

in Sh2-252 are comparable with those obtained for NGC 1893 (a = 0.34 ± 0.07), Stock 8 
(a = 0.31 ± 0.02) Be 59 (a = 0.27 ± 0.02) and NGC 1624 (a = 0.30 ± 0.06). 



3.6 Initial mass functions of sub-regions in Sh2-252 

One of the most fundamental disciplines of astrophysical research is the origin of stars and 
stellar masses. The distribution of stellar masses that form in a star formation event in a 
given volume of space is called IMF and is one of the most important measurable quantities 
in star formation studies. Together with star formation rate, the IMF dictates the evolution 
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and fate of star clusters and galaxies (Kroupa 2002). Young clusters and star forming regions 
are important objects to study the IMF since their mass functions (MFs) can be considered 
as IMFs as they are too young to loose a significant number of the members either by 
dynamical or stellar evolution. 

The MF is often expressed by the power law, A^(logm) oc mF and the slope of the MF 
is given as 

T = d log A^(log m)/d log m 

where A'^(logm) is the number of stars per unit logarithmic mass interval. For the mass 
range 0.4 < M/Mq ^ 10, the classical value derived by Salpeter (1955) is F = —1.35. 

One of the main objectives of this work is to measure the MFs of the sub-regions of 
Sh2-252. As discussed in section 13. 4. 4^ the majority of the sources lying above the 5 Myr 
isochrone of the CMDs could be probable PMS members of the region and most of them 
have a mass range of 0.3 - 2.5M0. We calculated the MFs of these PMS sources falling above 
5 Myr and mass range of 0.3 - 2.5M0 for the five sub-regions. To remove the contamination 
due to field stars from the PMS sample, we statistically subtracted the contribution of field 
stars from the observed CMDs of sub-regions (see Fig. [TTl) using the procedure described 
in our earlier studies (e.g., Pandey et al. 2008, Jose et al. 2008; 2011). After statistically 
subtracting the field star contribution, we used all the candidate PMS sources lying above 
the 5 Myr isochrone and having mass in the range of 0.3 - 2.5M0 to calculate the MFs. The 
MFs for the candidate PMS sources were obtained by counting the number of stars in various 
mass bins, shown as evolutionary tracks by Siess et al. (2000) in Fig. [HI Necessary corrections 
for data incompleteness as a function of magnitude were taken into account to calculate the 
MFs. The CFs given in Table [3] for regions C and NGC 2175s were used to correct for their 
data incompleteness whereas the average value of the CFs for regions C and NGC 2175s was 
applied to regions A, E and Teu 136 to correct for their data incompleteness. The resulting 
field star decontaminated and completeness corrected MFs of the PMS sources of the five 
sub-regions of Sh2-252 region are shown in Fig. [131 The slopes, F, of the MFs in the mass 
range 0.3 ^ M/Mq < 2.5 are found to be -1.33±0.51, -2.39±0.15, -1.40±0.30, -1.65±0.28, 
-1.82±0.26, respectively, for regions A, C, E, NGC 2175s and Teu 136. The slopes of the MFs 
of the regions A and E are in good agreement with that of the Salpeter (1955) value, and 
in view of large errors, those in the case of NGC 2175s and Teu 136 can also be considered 
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comparable to the Salpeter value. However, in the case of the region C, the MF seems to 
be steeper than the Salpeter value. In order to check this we also calculated the MF using 
J-band luminosity function (JLF). The procedure for this calculation has been mentioned 
in our earlier paper (Jose et al. 2011). The MF thus obtained for region C from JLF is 
shown in Fig. [13] as dashed line and the slope obtained is -2.10±0.54. Considering the large 
uncertainty in the IMF calculation, we conclude that all the regions have IMF slopes, which 
do not significantly deviate from the salpeter value. 

The shape of the stellar IMF and whether it is universal or not are key issues in as- 
trophysics. For clusters within 2 kpc, there is no compelling evidence for variations in the 
stellar IMF with respect to their environment (e.g. Meyer et al. 2000; Kroupa 2002; Chabrier 
2005). We have been pursuing studies of various star forming regions in the galaxy, hence a 
comparative study of their IMFs obtained using similar techniques will give useful informa- 
tion on IMFs. Our recent analysis of young clusters (age 2-4 Myr), viz., NGC 1893 (Sharma 
et al. 2007), Stock 8, NGC 1624 (Jose et al. 2008, 2011) and W5 E (Chauhan et al. 2011) 
have yielded values similar to the Salpeter value. In the present study, the MFs of all the 
sub-regions of Sh2-252 are found to be in general comparable to the Salpeter value. 

4 SUMMARY 

In this paper we studied the stellar contents of the H ll region Sh2-252 which itself contains 
four CH II regions such as A, B, C and E and two clusters NGC 2175s and Teu 136. We used 
the deep optical photometry in UBVRI bands, slit and slitless spectroscopic observations 
along with the JHK data from 2MASS for our analyses. An attempt has been made to 
identify and classify the massive members, to determine the fundamental parameters of 
Sh2-252 as well as to obtain the age and mass distribution of candidate PMS sources and 
finally to constrain the KLFs and IMFs of the sub-regions of Sh2-252. 

We have carried out optical spectroscopy of fifteen bright sources of the region, out 
of which, eight have been identified as massive members of spectral type earlier than B3. 
The CH II regions A, B, C and E are found to have at least one star of spectral type 
earlier than B3, whereas, the small optically bright cluster NGC 2175s is found to have 
four stars of spectral type earlier than B3 around its center. We have also identified the 
probable candidate ionizing sources of the CH ll regions A, B, C and E. From the spectro- 
photometric analysis, we derived the average distance of the region as 2.4 ± 0.2 kpc. The 
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Figure 13. MFs obtained from the optical data for the sub-regions of Sh2-252 after correcting for the field star contamination 
and data incompleteness. The </> represents Af/dlog m and the error bars represent ±\/N errors. The continuous lines show 
least-squares fits to the mass ranges described in the text. The value of the slopes obtained are given in each figure. The MF 
obtained from JLF for region C is shown in the second panel and the least-squares fit is shown by using the dashed line. 

distance estimated from the four bright stars of NGC 2175s is found to be in good agreement 
with the average distance of Sh2-252, which supports the association of NGC 2175s with 
the Sh2-252 complex. In this paper, we also report that the NIR embedded cluster Teu 136, 
located at the eastern edge of the complex as a sub-cluster of Sh2-252. Based on the optical 
colour-colour diagram and spectroscopic properties, the reddening E{B — V) of the massive 
members of the region is found to vary between 0.35 mag to 2.1 mag. 

Using slitless spectroscopy survey, we have identified 61 Ha emission sources in the re- 

© 2012 RAS, MNRAS 000, 000-000 



34 Jose et al. 

gion. We obtained optical photometry for 211 candidate YSOs of the region. The distribution 
of these candidate YSOs on the V/ {V — I) CMD shows that a majority of them have an age 
distribution between 0.1-5 Myr and masses in the range of 0.3 - 2.5 Mq. The stellar surface 
density distribution in K-hand shows that there is at least five sub-clusters, each associated 
with the regions A, C, E, NGC 2175s and Teu 136. The (J - H)/{H - K) C-C diagrams 
of the individual regions show that the regions A, C and E have more reddened population 
when compared to the clusters NGC 2175s and Teu 136, which are located towards the east 
of Sh2-252. The optical CMDs of the individual regions also show that the candidate PMS 
sources are distributed between 0.1-5 Myr, suggesting an age spread for them. 

We also calculated the KLFs of the sub-regions A, C, E, NGC 2175s and Teu 136 and 
have estimated as 0.31±0.03, 0.34±0.03, 0.25±0.05, 0.20±0.05 and 0.36±0.02, respectively. 
Within errors, the KLFs of all the sub-regions are similar and comparable to that of young 
clusters of age < 5 Myr. We also derived the MPs of PMS sample of the individual regions 
in the mass range of 0.3 — 2.5M0. The slopes of the MPs of all the sub-regions are found to 
match with the Salpeter value. 
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